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Die Aufgabe, welche iins beschaftigt, ist zu bestimmen, 
wie viel Ziffem bei einer jeden Anwe]idung des Newton- 
schen Verfabrens zuverlassig genau kommen werden. Ftir 
Gleichungen mit einer Unbekannten hat Fourier zur Beant- 
wortung dieser Frage foigenden Weg angegeben. 

Es sei die aufzulosende Gleichung 

fx = o. 
Die Wurzel a, welche man berechnen will, liege zwischen 
bekannten Grenzen a und b. Der Fall, dass a zugleich 
eine Wurzel der ersten oder zweiten Ableitung f x oder f 'x 
8ei, wird ausgeschloBsen und die Grenzen a und b mussen 
nahe genug sein, damit zwischen ihnen keine Wurzel von 
f X oder f"x liegt. Man hat 

f (a+a — a) = o . 
oder 

fa + (« — a) f' (a., a) = o, 
folglich 

— fa 
cc — a= 

f (a.. a); 

a.. a ist hier ein unbekannter zwischen a und a liearender 
Werth, eine Bezeichnung die in der Folge ebenso zu ver- 
stehen ist. Es ist also a — a ein spedeller Werth der Function 

— fa 
f X ' 
in welcher die Variable x von a bis b geht. Da man das 
Zeichen von f^x kennt, so kann man die extremen Werthe 
dieser Function unmittelbar berechnen. 

1* 



Es sei, abgesehen vom Zeichen, p der kleinste und q der 
grosste Werth der Function 

— _fe 

f X 

innerhalb der Grenzen a und b der Variablen x; ebenso 

sei r der kleinste und s der grosste Werth der Function 

— jfb 
fx~' 
Man hat dann folgende vier neue Paarė von Grenzen 

a + p und b + r, 
b + s und a + q, 
a + p und a + q, 
b + s und b + r. 
Fourier hat nun das Intervall der neuen Grenzen ermittelt 
und verglichen mit dem anfanglichen Intervall b — a. 

Das Intervall der Grenzen a + p und b + r erhalt den 
Ausdruck 

f" (a. .b),, ,2 

wenn 

(f a)2 > (i'hy, 
denn andernfalls vvare der Divisor f b. Fiir das Intervall 
von b -f- s und a + q besteht der namliche Ausdruck , nur 
ist der Divisor f 'b, wenn jener fa und fa, wenn jener 
f b ist. 

Die Grenzen a + p und a + q geben ein Intervall 

f' (a.. b) ,, ,. 

[For:7b)p ^^-*>^*- 

Fourier verlangt, lun die extremen Wertlie von f x be- 
rechnen zu konnen, dass auch die dritte Ableitung i***x in 
dem Intervall a.. b keine Wurzel habe. Bezeichnet dann 
f den numerisch grossten Werth der zweiten Ableitung 



sowie f den numerisch kleinsten Werth der ersten Ableitung 
innerhalb der Grenzen a und b der Variablen, so wird 

f" 



2f' 

numerisch grosser sein als das Intervall der Grenzen a + p 
und b + r. Fourier wendet nun das Naherungsverfahren 
auf diejenige der beiden anfanglichen Grenzen an, fiir welche 
die erste Ableitung numerisch am grossten ist ; er berechnet 
also, wenn dies die Grenze a, 

fa 

und bestimmt mit Hiilfe des Werthes 



2f' 

wie viel Ziffem bei seiner Operation zuverlassig genau 
kommen werden. Die letzte genaue Ziffer vermehrt oder 
vermindert er um eine Einheit und hat somit neue engere 
Grenzen a, und b,. Die Genauigkeit von 

fa, 



a. 



fa, 
wird wieder durch 

bestimmt, so dass der Factor 

f' 
2f 

den ganzen Verlauf der Rechnung regulirt. 

Man kann dieses Verfahren von Fourier auf Gleichungen 

mit zwei oder mehr Unbekannten ausdehnen; ich habe dies 

in meiner Inauguraldissertation gethan. Mehr jedoch als 

eine solche Ausdehnung scheint sich ein Weg zu empfehlen, 

den ich hier angeben will. 



Um zunachst eine Gleicbung zu betrachten, so sei wiedeT 

fx = o 
die aufzulosende Gleichung, a die zu berechnende Wurzel 
undaein bekan]iterNaherung8werth derselben. Man berechne 
h aus der Gleichung 

fa 4- bf'a == o. 
Es ist 

f (a + h + a — a — h) = o 
oder 

f(a + h) + (a — a — h)f' (a+h..a)=o, 
also 



a — a — h = 



— f(a + h) 



f'(a + h..a)' 
und da 

f(a+h) = ^-:^i^±^U^ 

80 folgt 

2r(a + h..a) 
Dieser Ausdruck wird uns dažu dienen, um den Grad der 
Genauigkeit von a + b zu bestimmen. 

Es moge die Wurzel a eingeschlossen sein in zwei 
Grenzen b und c. Indem w den Fall von fa = o aus- 
schliessen, verlangen wir diese so engė, dass die Funetion i* 
innerhalb dieser Grenzen der Variablen nicht yerschwindet. 



Femer sei M der numeriscfa grosste Werth, welchen die 
Function 

— f^y 
2fx 

anniinmt, wahrend die Variablen x und j zwischen b und c 
bleiben. Wenn niin b und c einander so nahe sind, dass 

M (b— c)2 

Yon einer hoheren Decimalordnung ist als b — c, so wird 
um 80 mehr 

Ton einer hoheren Ordnung sein als h, sobald man fiir h 
einen numerisch kleinem Werth setzt als b — c. 

Wenn wir also dann mit einem Werthe a zwischen 
b und c das Newton'sche Verfahren machen d. h. 

berechnen, und wenn der so erhaltene neue Werth gleich- 
falls zwischen b und c liegt, so muss er nothwendig eine 
grossere Anzahl genauer Decimalstellen haben als a. Denn 
er wurde zur volligen Genauigkeit erst in einer hoheren 
Stelle einer Correction bediirfen als diejenige ist, mit welcher 
h anfangt. Gesetzt nun es lage a + h ausserhalb des Inter- 
yalls b.. c aber unendlich nahe bei einer der Grenzen, so 
wurde a -f- h noch immer genauer sein miissen als a. Da- 
raus folgt dass wenn man fiir a die der Wurzel nahere 
Grenze nimmt, dann gewiss a+h innerhalb der Grenzen 
liegen wird. Auch die entfemtere Gbrenze kann uns nicht 
aus b.. c hinaus fiihren, wenn 

M(b— c)2 

kleiner ist als das Intervall zwischen a und der n'achsten 
Grenze. 
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Sobald niin 

M(b — c)2 

von einer hoheren Decimalordnung ist als b — c, so kann 
man zur Approximatioii schreiten und wird eine stets wach- 
sende Anzahl genauer Decimalziffem erhalten. AUerdings 
kann es geschehen, wenn namlich die Wurzel der einen 
Grenze sehr viel naher liegt als der anderen und man dies 
nicht erkannt hat, dass indem man die Naherung mit der 
entfemteren Grenze macht, man iiber die andere hinaus 
kommt. Aber dann wahle man diese andere Grenze, so 
kann dies nicht geschehen. 

Zur Beurtheilung der Approximation ist es zunachst 
nothig, dass man die beiden extremen Werthe der Grosse 

— ^Pį 
2f'x 

in dem IntervalI b . . c fiir x und y kenne. Sie seien M und 
N und zwar M der algebraisch grossere, N der algebraisch 
kleinere. Es besteht dann die Ungleichung 

Nh2<a — a— h<Mh2. 

Will man also mit einem Werthe a zwischen b und c die 
Approximation machen, so berechne man von 

h = -į - 
i'a 

yorlaufig nur die erste Decimalziffėr, welche hi und ver- 

grossert um eine Einheit h2 heissen mag. Aus den vier 

Zahlen 

M N hi h2 

bilde man nun die beiden Correctionen fiir 

a — a — h 
und setze die Berechnung von 

h = -^ — 
ra 



nicht weiter fort als bis zu der Stelle, welche durch diese 
Correctionen alterirt wirdy wodurch man unmittelbar die 
beiden neuen und engeren Grenzen der Wurzel erhalt. 

Die hier angegebene BestimmuDg der wirklichen An- 
naherung ist von der Fourier'schen dadurch verechieden, 
dass Fourier mit dem genaherten Werthe a + b oder 

fa 

die andere Grenze 

b — ^- 
ra 

verbindet, wogegen wir die extremen Werthe der Function 

— f(a + h) 
fx 

berechnet und dadurch Grenzen fiir 

a — a — h 
erhalten haben. Als Mass der Approximation ergibt sich hier 

2fx 



und bei Fourier 






2f 

£s musste oben vorausgesetzt werden, mit Ausschliessung 
von f a = o , dass die erste Ableitung f'x keine Wurzel 
innerhalb der anfanglichen Grenzen habe. Ob die zweite 
Ableitung f" eine Wurzel in diesem Intervall hat, ist da- 
gegen gleichgultig, da es nur auf die numerische Klein- 
heit von 

L h2 

2f' 

ankommt. Will man das Sturm'sche Verfahren nicht auf 
fx anwenden9 so kann man sich noch auf folgende Weise 
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iiberzeugen, dase f'x nicht versch^indet. Die Grenzen 
b iind c, zwischen welchen x zu nehmen ist, werdeii das 
namliche Vorzeichen haben. Indem man also die positiven 
und die negativen Glieder besonders zusammenfasst, wird 

f'x = qpx — tfjx. 
Wenn nun 

qpb — ifjc 
und 

q)C — i^b 
einerlei Vorzeichen haben, so ver8chwindet sicher f'x nicht. 
Und wenn {'a nicht Nuli ist, so kann man b und c immer 
so nahe zusammenziehen, damit beide 6renzwerthe von 
cf — tĮj einerlei Zeichen erhalten. Bei Functionen mit 
mehreren Variablen muss man ohnehin dieses Mittel an- 
wenden, wenn man sich iiberzeugen will, dass eine solche 
Function innerhalb gewi88er Grenzen ihrer Variablen nicht 
ver8chwindet. Aus 

(fh yfh qp'b 1^'b 

(po tpc C/'C 1^'c . 
kann man nun die extremen Werthe M und N von 

unmittelbar berechnen, welche die Annaherung an die Wur- 
zel reguliren vermoge der obigen Ungleichung 

Nh2<a — a— h<Mh2. 



Wir gehen nun iiber zu dem Faile zweier simultanen 
Gleichungen. Diese seien 

f =s o, C) = o 

und die zu berechnenden Wurzel]i a und /9, so dads also 

iaP = Oy qpa/?=o. 
Sind a und b bekannte Nahertuigswerthe, so ziehe man 
Annaherungen h und k aus den beiden linearen Gleichungen 

fab + hf ab + kf,ab = o 
9)ab + hc/'ab + kijp^ab = o, 
woduFch man die neuen Werthe 

a -{- h == a, und b + k == b, 
erhalt. Wir suchen nun wieder 

rt — a, und /? — b, 
in Grenzen einzuschliessen, zu welchem Zwecke die Glei- 
chungen dienen 

f(a, + a — a,, b, + /? — b,) = o 
qp(a, + a — a,, b, + /? — b,) = o, 
aus welchen folgt 

o = fa,b, + (a — a,) i* (vL,..a, b,../9) 

+ (/9— b,)f,(a,..a,b,../?) 

o = (ps.,h, + (a — a,) (p* (eL,..ay h,. ./J) 

+ O*— M V; (a, . . a, b, . . /?); 

hier sind a,.. a und b,../9 unbekannte zwi8chen den be- 

zeichneten Grenzen liegende Werthe, aber nicht dieselben 

in den Ableitungen von f wie in denen von y. Setzen wir 

a, . . a = X und b, ../} = j 
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in den Ableitungen von f, dagegen 

a, . . a = I und b, . . /S = iy 
in den Ableitungen von qp. Aus den vorstehenden (jlei- 
chungen ergeben sich dann 

Es bleiben nun noch auszudriicken fa, b, und c)a,b,. F^*^ 
diese hat man 

fa,b, = 2^ f x,y, + hkf;x,y, + j i„x,y, 

h2 k^ 

qpa;b, = ^ qp"|,i?, + hk^;!,??, + y (p„į,v,l 

hier sind x, und |, so wie y, und i;, neue unbekannte Werthe, 
jene beiden zwi8chen a und a 4~ hj diese beiden zwi8chen 
b und b + k. Setzt man die Werthe fiir f a, b, und qpa, b, 
ein in a — a, und /? — b,, so erhalt man 

a-a— h = Ah2 + Bhk + Ck^ 
/9— b— k = 8lh2 + ©hk + ek2. 

In den Coefficienten sind enthalten die acht Grossen x und I9 
welche liegen zwischen a + h und a, y und tį zwi8chen 
b + k und /?, X, und į, zwischen a und a + b, y, und i], 
zwischen b und b + k. 

Es sei nun die Wurzel a eingeschlossen in den be- 
kannten Grenzen c und d, die Wurzel fi in den Grenzen 
e und g. Den Fall wo 

t*a/}(p,a/3 — f,afi(p'afi = o 
schliesse ich aus und verlange die Grenzen c und d, e und 
g so engė, dass der gemeinschaftliche Nenner 

iunerhalb der Grenzen c und d fur x und |, e und g fur 
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y und fj nicht verschwmde, wa8 durch das friiher ange- 
deutete Mittel zu beurtheilen ist. 

Man berechne nun die extremen Werthe der sechs 
Coefficienten A.. 8.. einmal durch Vergrosserung der po- 
sitiven und Verkleinerung der negativen Glieder und dann 
umgekehrt durch Verkleinerung der positiven und Ver- 
grosserung der negativen Glieder, wobei den Grossen xįx,į, 
die Werthe c und d beizulegen sind, den Grossen yriytVf 
die Werthe e und g. 

Die Buchstaben A*. S(.. mogen jetzt die numerisch 
grossten Werthe dieser Coefficienten bezeichnen positiv ge- 
nommen. Die unbekannte Entfemung der Wurzel a von 
ihrer nachsten Grenze sei 8 und der Wurzel /? von der 
ihrigen €. Wenn dann die beiden Ungleichungen erfullt sind 

Ah2 + Bhk+Ck2<3 
ah2 + S3hk + Sk2<6, 

wo h =r c — d und k = e — g positiv zu nehmen ist, so 
wird jeder Punkt ab in dem Bechteck c e., d g durch das 
Newton'sche Verfahren zu emem Punkte a + 1^9 b + k 
fiihren, welcher gleichfalls in demselben liegt. 

Bezeichnen m und n die links stehenden Werthe in den 
beiden Ungleichheiten , so werden wir ihrer Erfullung sicher 
sein, wenn die Grossen 

c + m und d — m 

e + n und g — n 
noch Grenzen sind fur a und /?• Wenn m von hoherer 
Decimalordnung ist als c — d und n von hoherer als e — g, 
so muss wenigstens eine Wurzel viel naher an einer von 
ihren beiden Grenzen liegen als an der anderen, wenn 

c + m und d — m 

e + n und g — n 
keine Grenzen mehr sein soUten. 
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Um jedoch auch diesen Fall einzuschliessen, so nehme 
man die erweiterten Grrenzen 

c — fA und d~\- fi 
e — v und g + v, 
wo fz und v je eine Einheit sein mogen in den respectiven 
SchlusBstellen von c und d, e und g. Man priife nun, ob 
fur diese erweiterten Grenzen die Werthe m und n respec- 
tive kleiner ausfallen als ^ und i/. Wenn dies der Fall, 
80 besitzt man in den erweiterten Grenzen ein Sechteck, 
aus welchem keine Entfemung bei beliebiger WiederhoIung 
des Newton'8chen Verfahrens mogUch ist. 

Wenn /u > v, so wird bei einer hinreichenden Kiein- 
beit von c — d zwar m <C fi} ob abcr auch n < v hangt 
davon ab, dass e — g nicht zu viel kleiner sei als c — d. 
Allein, wenn man die eine Wurzel in sehr viel engere 
Grenzen eingeschlossen hat als die andere, so kann man 
sich durch Erweiterung der zu engen Grenzen helfen, ohne 
dass man darum die Genauigkeit aufgiebt, die man schon 
hat, wie aus dem Folgenden zu ersehen ist. 

Wir mūssen nun voraussetzen, dass man sich Grenzen 
verschafit habe, reiche den beiden oben stehenden Un- 
gleichungen geniigen, was auf die angegebene Art zu er- 
kennen ist. Bezeichnen wir jetzt die algebraisch grossten 
Werthe der sechs Coefficienten mit A.. 31.., die algebraisch 
kleinsten mit Ao..$(o.*- Es bestehen dann fiir die Beur- 
theilung der Approximatton die Ungleichungen 

, > Aoh2 + Bohk + Cok2 
a_a— h ^ ^į2 ^ Btk ^ ck2 

^ , , >8loh2 + ©ohk-f 6ok2 
^_b — k <;8ih2 + a3hk+6k2 

bei gleichen Zeichen von h und k; bei ungleichen Zeichen 
vertauschen sich Bo und B, JBo und S3. 
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Vergleicht man diese mit den beiden obigen, so zeigt 
sich dass wenn numerisch 

k = e — g 

dann nothwendig a — a — h von hoherer Decimalordnung 
ausfallt als h; wogegen /? — b — k von hoherer Ordnung 
wird als k y wenn 

h ^ c— d 

k < e— g 

Eine wirkliche Annaherung iSndet daher nothwendig stati 

wenig8tens mit einer der beiden Grossen h und k. Sie 

kann freilich auch mit beiden stattfinden. 

Um nun die Genaoigkeit von a + h und b + k zu be- 
Btimmen, so berechne man von h vorlaufig nur die erste 
Decimalstelle 9 die hi und um eine Einheit vergrossert ht 
heissen mag, und ebenso ki und k2. Mittelst der obigen 
Ungleiohheiten bilde man jetzt die Grenzen fur a — a— h 
und /? — b — k. Fallen beide Grenzenpaare von hoherer 
Ordnung aus als respective h und k, so berechne man beide 
Werthe a + h und b + k bis zu den Stellen, welche durch 
ihre Correctionen alterirt werden , andemfalls nur den einen 
Werth dessen grossere Genauigkeit verbiirgt ist. Man er- 
halt dann mittelst der Correctionen engere Grenzen fur eine 
Wurzel oder fiir beide. Es ist jedoch leicht zu seben, 
dass die Annaherung nicht immerfort an dieselbe Wurzel 
stattfinden kanu, ohne dass sich eine Annaherung an die 
andere Wurzel er^be. Es ist nicht einmal moglich sich 
der einen Wurzel immerfort zu nahem, wahrend man die 
Annaherung an die andere absichtlich unterlasst. Man wd 
somit die Annaherung an beide Wurzeln bis auf jeden belie- 
bigen Grad treiben konnen und eine immer wachsende 
Anzahl genauer DecimalziiFern erhalten. 
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Die Erstreckung unseres Verfahrens auf Gleichungen 
mit mehr als zwei Unbekaimten ergibt sich jetzt yon selbst ; 
mk einer von den Grofisen hkl.. wirc[ nothwendig eine 
Ąnnahening stattfinden, wie man sich leicht uberzeugt, 
und es kaan keine Wurzel bestandig von der Annaherung 
ausgeschlossen bleiben. Auch sieht man, in wie weit es 
bei transcendenten Gleichungen Anwendung finden kaim. 



Entes BeispieL 

£s sei 2u berechnen diejenige Wiirzel der Gleichung 

x»— 2x — 5 = 0, 

welche liegt zwischen 2 und 2,1. 

Der Regulator wird fiir dieses Intervall 

— 0,63h2 < a — a — h < --0,53h2. 

Wenn man hier die beiden Grenzen um M (b — c)^ = 0,0068 
verengen wollte, so wurde man die Wurzel nicht mehr 
zwischen ihnen haben. Man ist daher nicht versichert in 
dem Intervalle zu bleiben. In der That liefert die Ann'ahe- 
rung mit der Grenze 2 gerade die andere Grenze 2,1, 
Fiir a = 2,1 erhalt man nun 

— 0,006 < h < — 0,005, 
folglich 

— 0,000023 < a — a — h < —0,00001. 

Nun wird 

a -+- h = 2,094568 ... , 

also sind die neuen Grenzen 

2,09454 und 2,09456. 

Bei der zweiten Approximation wird fiir a = 2,09454 

0,000011 < h < 0,000012, 
folglich 

—0,01091 < a — a— h < -0,Oio64. 
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Da sich nun ergibt 

a + h = 2,09455 14816 16.., 
8o werden die neuen Gr^izen 

2,09455 14815 2 
2,09455 14815 6. 



Zweite8 BeispieL 

Es seien aufzulosen die beiden GleichuDgen 

x3 — y3 — 2x* + 4xy = o 

xy2 -f-x2 — y2 = o, 

und zwar soll man berechnen die beiden Auflosungen a 
und fi, welche liegen respectire zwischen 

0,77 und 0,78 

1,62 und 1,63. 

Der Regulator fiir a + b ist in diesem Intervall 
(—0,33.. —0,24) h2 + (-0,88.. —0,69) hk 
+ (-0,20..— 0,11) k2, 

der Regulator fiir b + k 

(— 0,28 . . — 0,20) h2 + (— 0,04 . . + 0,02) h k 
+ (—1,31..— 1,06) kl 

Wenn man bier die fruher durch m und n bezeichneten 
Werthe berechnet, welche 0,000141 und 0,000163 werden, 
und die Grenzen um eie verengt, so behalt man noch 
Grenzen, so dass man sicher ist aus dem anfanglichen 
Intervall nicbt herauszukommen und der Regulator immer 
derselbe bleibt. 

Bei der ersten Approximation wird nun fiir 

a = 0,77 und b = 1,62 
0,003 < h < 0,004 
0,005 < k < 0,006, 
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daher 

—0,00004 < a -a, < —0,00001 
— 0,00006 < /? - b, < — 0,00002. 

Da sich nun ergibt 

a, = 0,77859 . . 

b, = 1,62572 . . , 
80 sind die neuen Grenzea fiir a 

0,77355 und 0,77359, 
fur /? 

1,62566 und 1,62571. 
Bei der zweiteD Approsimadon wird, indem man von 
den kleinsten Werthen ausgefat, 

0,000021 < į< 0,000022, 

somit 

— 0,0''7 < a - a, < — 0,0«4 

-0,0»8 < /9— b, < — 0,0«5. 
Nun ergeben sich 

a, = 0,77357 17770 . . 

b, = 1,62568 10257.., 
daher werden die neuen Grenzen fiir a 

0,77357 17763 
0,77357 17767 



und fiir /? 



1,62568 10249 
1,62568 10253. 
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